Field sweep rate dependence of the coercive field of single-molecule magnets: a 
classical approach with applications to the quantum regime 
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A method, based on the Neel-Brown model of thermally activated magnetization reversal of a 
magnetic single-domain particle, is proposed to study the field sweep rate dependence of the coercive 
field of single-molecule magnets (SMMs). The application to Mni2 and Mns4 SMMs allows the 
determination of the important parameters that characterize the magnetic properties: the energy 
barrier, the magnetic anisotropy constant, the spin, to, and the crossover temperature from the 
classical to the quantum regime. The method may be particularly valuable for large SMMs that do 
not show quantum tunneling steps in the hysteresis loops. 

PACS numbers: 75.50.Xx, 75.75.+a, 75.45.+j, 75.50.Tt 



I. INTRODUCTION 

Single-molecule magnets (SMMs) exhibit slow magne- 
tization relaxation of their spin ground state, which is 
split by axial zero-field splittingi* 2 ^*^. They are now 
among the most promising candidates for observing the 
limits between classical and quantum physics since they 
have a well defined structure, spin ground state and 
magnetic anisotropy. An important effort in synthetic 
chemistry has led to a quickly growing number of SMMs 
with an increasing number of magnetic centers. Recently, 
the molecular (or bottom-up) approach has reached the 
size regime of the classical (or top-down) approach to 
nanoscale magnetic materials-®. Indeed, a giant Mn 8 4 
SMM was reported with a 4 nm diameter torus struc- 
ture, exhibiting both magnetization hysteresis and quan- 
tum tunneling. 

The study of such large systems is greatly complicated 
by the fact that the spin Hilbert space is huge and it is 
impossible to treat such systems with exact matrix diago- 
nalization methods. However, since some SMMs are now 
as large as some classical nanoparticles, it raises the in- 
teresting possibility that classical models commonly em- 
ployed to study the latter may be used to obtain a first- 
order understanding for large molecular systems. Indeed, 
we herein propose and demonstrate the use of the classi- 
cal Neel-Brown modeliS*^ of thermally activated magne- 
tization reversal of a magnetic single-domain particle in 
order to study large SMMs. The proposed method allows 
us to determine important parameters that characterize 
the magnetic properties of the SMM: the energy bar- 
rier, the magnetic anisotropy constant, the spin, to, and 
the crossover temperature from the classical to the quan- 
tum regime. The method is particularly useful for SMMs 
having low-lying energy states and not showing quantum 
tunneling steps in hysteresis loops. In such systems elec- 
tron paramagnetic resonance (EPR) measurements often 
exhibit only very broad absorption peaks which do not 
allow the determination of the magnetic anisotropy. 

In this letter, we apply the method to two systems 
whose properties are already known: (i) a Mni2 SMM 



with a well-characterized S = 10 ground state and (ii) a 
giant Mn§4 SMM with S — 6 which has many low-lying 
energy states with higher spin values. 



II. NEEL-BROWN MODEL 

The presented method is based on the Neel-Brown 
model of thermally activated magnetization reversal of 
a magnetic single-domain particle which has two equiv- 
alent ground states of opposite magnetization separated 
by an energy barrier due to magnetic anisotropjiSiS. The 
system can escape from one state to the other either by 
thermal activation over the barrier at high temperatures 
or by quantum tunneling at low temperatures. At suf- 
ficiently low temperatures and at zero field, the energy 
barrier between the two states of opposite magnetization 
is much too high to observe an escape process. However, 
the barrier can be lowered by applying a magnetic field 
in the opposite direction to that of the particle's magne- 
tization. When the applied field is close enough to the 
reversal field, thermal fluctuations are sufficient to allow 
the system to overcome the barrier, and the magnetiza- 
tion is reversed. 

This stochastic escape process can be studied via the 
relaxation time method consisting of the measurement of 
the probability that the magnetization has not reversed 
after a certain time. In the case of an assembly of iden- 
tical and isolated particles, it corresponds to measure- 
ments of the relaxation of magnetization. According to 
the Neel-Brown model, the probability that the magne- 
tization has not reversed after a time t is given by: 



P(t) = e 



- „-t/r 



(1) 



and t (inverse of the reversal rate) can be expressed by 
an Arrhenius law of the form: 



T(T,H) = T e AE W/ k * T 



(2) 



where AE(H) is the field dependent energy barrier height 
and tq is the inverse of the attempt frequency. In most 
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FIG. 1: (color online) Coercive field H c for Mni2 as a function 
of (a) temperature and (b) field sweep rate. Note the steps of 
H c coming from the resonant tunneling steps in the hysteresis 
loops. 



cases, AE(H) can be approximated by 
AE(H) =E Q (1- H/H° c ) 



(3) 



where H® is the reversal field at zero temperature, Eq is 
the barrier height at zero applied field, and a is a con- 
stant of the order of unity (for most cases 1.5 < a < 2). 
In the case of a Stoner-Wohlfarth particl e 10 ! 11 with uni- 
axial anisotropy and the field applied along the easy axis 
of magnetization, all constants can be determined an- 
alyticalljii .: a = 2, E = KV, and # c ° = 2K/M S , 
where K is the uniaxial anisotropy constant, V is the 
particle volume, and M s is the saturation magnetization. 
For SMMs with dominating uniaxial anisotropy: a — 2, 
Eq = DS 2 , and H® = 2DS/ g^iQ^B- However, in gen- 
eral, all constants depend slightly on fine details of the 
magnetic anisotropy and the direction of the applied field 

In order to study the field dependence of the relax- 
ation time t(T, H) and to obtain the parameters of the 
model, the decay of magnetization has to be studied at 
many applied fields H and temperatures T. This is ex- 
perimentally very time consuming and complicated by 
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FIG. 2: (color online) Coercive field H c for Mn84 as a function 
of (a) temperature and (b) field sweep rate. 



the fact that the equilibrium magnetization is tempera- 
ture dependent and difficult to obtain for long relaxation 
times. In addition, for fast relaxation times the initial 
magnetization depends on the field sweep rates to apply 
the field. The number of exploitable decades for r values 
is therefore limited for relaxation time measurements. 

A more convenient method for studying the magneti- 
zation decay is by ramping the applied field at a given 
ratei^ and measuring the coercive field H c (the field value 
to obtain zero magnetization) , which is then measured as 
a function of the field sweep rate and temperature. 

The mathematical transformation from a reversal time 
probability (Eqs. Hand[2J to a reversal field probability 
was first given by Kurkijarvii^ for the critical current in 
SQUIDs. Later, Gunther and Barbara calculated simi- 
lar expressions for magnetic domain wall junctions^. A 
more general calculation was evaluated by Gargit Here, 
we use a simplified version (see annex)i4 and approxi- 
mate the mean reversal field of an assembly of identical 
particles or SMMs by the coercive field H c : 



H c (T,v)^H° c [l 



kT . 



■In 



(;) 



(4) 



3 



where the field sweeping rate is given by v = dH/dt; 
H® is the coercive field at zero temperature, and c 
depends on the details of the approximations: c = 
H°k B T/[T aE {l - HJH^- 1 ] in referenced, c' = 
(Eo/kT) 1 / /(root) in reference^, and it can be taken 
constant when the exact value of tq is not needed. We 
applied the three approximations to nanoparticlesia and 
here to SMMs and found that the first approximation 
gives a tq which is closest to that extracted from an Ar- 
rhcnius plot. 



III. APPLICATION TO MN t2 AND MN 84 

In the present work, the method 
is applied to two SMMs: (i) 
[Mni 2 Oi 2 (0 2 CCH 2 Bu%(CH 3 OH) 4 ] • CH 3 OH and 

(ii) [Mn 84 72 (0 2 CMe) 78 (OMe) 2 4(MeOH) 12 (H 2 0)42(OH) f 
called respectively Mni 2 and Mn§4 henceforth. Full 
details of the synthesis, crystal structure and magnetic 
characterization are presented elsewher o 6 ! 18 . 

The magnetization measurements were performed by 
using (i) a magnetometer consisting of several 6x6 fim 2 
Hall-bars^ and (ii) an array of micro-SQUIDs on top 
of which single crystals of Mni 2 and Mn§4 were placed, 
respectively. The field can be applied in any direction 
by separately driving three orthogonal superconducting 
coils. The field was aligned with the easy axis of magne- 
tization using the transverse field method^. 

Typical hysteresis loops of both systems can be found 
elsewhere. Mni 2 displays hysteresis loops with a series 
of quantum steps separated by plateaus^i whereas Mn§4 
shows a smooth hysteresis without steps^. In order to ap- 
ply the above method, the temperature and field sweep 
rate dependences of the coercive fields H c were measured 
and plotted in Figs ^ an d [21 Note that the curves for 
Mni 2 show steps coming from the steps in the hysteresis 
loopsSi. As expected for a thermally activated process, 
H c increases with decreasing temperature and increas- 
ing field sweep rate. Furthermore, all our measurements 
showed an almost logarithmic dependence of H c on the 
field sweep rate (Figs.^ 3 and[2t>). H c becomes tempera- 
ture independent below about 0.6 and 0.3 K, respectively 
for Mni 2 and Mn§4. 

The validity of Eq. 0] was tested by plotting the set 
of H c (T,v) values as a function of [Tln(cfv)] 1 ' 2 where 
c = H°k B T/T 2E (l - H c /H°). If the underlying model 
is sufficient, all points should collapse onto one straight 
line by choosing the proper values for the constant To. 
We found that the data of H c (T, v) fell on a master curve 
provided To = 2.1 x 10~ 7 s in Fig. Inland 2 x 10~ 7 s in 
Fig.0] Whereas for Mn§4 the master curve is straight, it 
presents steps for Mni 2 . 

At low temperatures, strong deviation from the master 
curves are observed. In order to investigate the possibil- 
ity that these low-temperature deviations are due to es- 
cape from the metastable potential well by tunneling, a 
common method for classical models is to replace the real 
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FIG. 3: (color online) (a) Scaling plot of the coercive field 
H C (T, v) of Mni2 for field sweep rates between 0.0001 and 0.1 
T/s and several temperatures: 0.1 K, 0.4 K, from 0.6 to 1 
K in steps of 0.05 K, and from 1 to 4 K in steps of 0.1 K. 
The arrows indicate the step index n = — (m + m') where m 
and m' are the quantum numbers of the corresponding level 
crossing. Note the parity effect of the steps: even n have 
larger steps than odd n. (b) Same data of H C (T, v) and same 
scales but the real temperature T is replaced by an effective 
temperature T* (see inset) which restores the scaling below 
1.1 K. 



temperature T by an effective temperature T*(T) in or- 
der to restore the scaling plot 22 . In the case of tunneling, 
T*(T) should saturate at low temperatures. Indeed, the 
ansatz of T*(T) as shown in the inset of Figs. 12b and^Jj, 
can restore unequivocally the scaling plot demonstrated 
by a straight master curve (Figs. [2t> and^>). The flat- 
tening of T* corresponds to a saturation of the escape 
rate, which is a necessary signature of tunneling. The 
crossover temperature T c can be defined as the temper- 
ature where the quantum rate equals the thermal one. 
The inset of Figs. 03 and Ha gives T c = 0.97 and 0.37 
K for Mni 2 and Mn§4, respectively. The slopes and the 
intercepts of the master curves give E = 72.4 and 15.6 
K and = 9.1 and 3.1 T, respectively for Mni 2 and 
Mn§4. The Eq values are in good agreement with those 
extracted from Arrhenius plots (69 K and 18 K for Mni 2 
and Mn§4, respectively)^^. This result allows us to esti- 
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FIG. 4: (color online) (a) Scaling plot of the coercive field 
H c (T,v) of Mn84 for field sweep rates between 0.001 and 
0.1 T/s and several temperatures (Fig. |5] (b) Same data 
of H c (T,v) but the real temperature T is replaced by an ef- 
fective temperature T* (see inset) which restores the scaling 
below 0.5 K. 



mate the spin ground state using S = 2E /(g^B^oHc )■ S 
— 11 and 7 for Mni2 and Mn§4, respectively. This differs 
slightly from S — 10 and 6 determined via magnetization 
measurements. This deviation is due to quantum effects 
in the thermally activated regime and is considered fur- 
ther below. 

Several points should be mentioned: (i) the classi- 
cal regime of the model corresponds in most SMMs to 
the thermally activated tunneling regime with tunneling 
close to the top of the energy barrier. Because all param- 
eters are deduced from this regime, small deviations from 
the exact values are expected; (ii) Eq. 0] is not valid for 
fields which are close to H = because the model only 
takes into account the transitions from the metastable 
to the stable well. However, close to H = 0, transitions 
between both wells are possible leading to a rounding of 
the master curve at small fields; (in) the method can be 
applied to powder samples with random orientations of 
the molecules. In this case, a « 1.5, vEq — DS 2 where 
v can be calculatedi2*i£, and the intercept of the mas- 
ter curve gives H®/2; (iv) in the case of a distribution 



FIG. 5: (color online) (a) Field dependence of the energy 
barrier of Mni2 obtained from Eq. [H]and the set of H c (T,v) 
data from Fig. The arrows indicate the step index n = 
— (m + m') where m and m' are the quantum numbers of the 
corresponding level crossing. Note the step like reduction of 
the energy barrier due to resonant tunneling and the parity 
effect of the steps: even n have larger steps than odd n. The 
dotted line gives the classical barrier AE — Eo(l — H/H&) 2 
with E = 74 K and H a = 9.8 T. (b) Same data of H C (T, v) but 
the real temperature T is replaced by an effective temperature 
T* (see inset of Fig-Efc). 



of anisotropics, different parts of the distribution can be 
probed by applying the method at different M values; 
(v) this method is insensitive to small intermolecular in- 
teractions when H c is larger than the typical interaction 
field; and (vi) the method can be generalized for ID, 2D, 
and 3D networks of spins. In this case, Eq. [31 describes a 
nucleation barrier. 



IV. CONCLUSION 

We introduce a method that uses the temperature and 
field sweep rate dependences of the coercive field of large 
SMMs in order to determine several parameters that 
characterize the magnetic anisotropy and the relaxation 
dynamics. We have succesfully applied the method to 
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FIG. 6: (color online) (a) Field dependence of the energy 
barrier of Mns4 obtained from Eq. [S]and the set of H c (T,v) 
data from Fig. 0] (b) Same data of H c (T,v) but the real 
temperature T is replaced by an effective temperature T* (see 
inset of Fig. gjj). 



two test cases: (i) a Mni2 SMM with well known mag- 



netic properties; and (ii) a giant Mn§4 SMM which has 
many low-lying energy states with higher spin values. We 
believe that this method is an important tool to charac- 
terize magnetically new molecular systems of great com- 
plexity which do not allow a detailed undestanding on 
the quantum level. 

This work was supported by the EC-TMR Network 
QuEMolNa (MRTN-CT-2003-504880), CNRS, Rhone- 
Alpes funding, and NSF. 



V. ANNEX 

The probability density of reversal of a stochastic pro- 
cess is 

dt t . w 

and the maximum of the probability density can be de- 
rived from 



d 2 P P f l+ dr 



dt 2 



dt 



(6) 



This gives the following general result for the maximum 
of the probability density 



dt 

The application of the result to Eq. [21 leads to 



AE(H) = fc B Tln 



dE dH 
, T °dH dt . 



(7) 



(8) 



Using Eqs.EJand|U we find Eq.H 

Eq.[S]can be used to plot directly the field dependence 
of the energy barrier (Figs. |SJ and |SJ 



1 R. Sessoli, H.-L. Tsai, A. R. Schake, S. Wang, J. B. Vin- 
cent, K. Folting, D. Gatteschi, G. Christou, and D. N. 
Hendrickson, J. Am. Chem. Soc. 115, 1804 (1993). 

2 R. Sessoli, D. Gatteschi, A. Caneschi, and M. A. Novak, 
Nature 365, 141 (1993). 

3 G. Christou, Polyhedron XX, XXX (2005). 

4 N. Aliaga-Alcade, R.S. Edwards, S.O. Hill, W. Wernsdor- 
fer, K. Folting and G. Christou, J. Am. Chem. Soc. 126, 
12503 (2004). 

5 N. E. Chakov, M. Soler, W. Wernsdorfer, K.A. Abboud, 
and G. Christou, Inorg. Chem. 44, 5304 (2005). 

6 A.J. Tasiopoulos, A. Vinslava, W. Wernsdorfer, K.A. Ab- 
boud, and G. Christou, Angew. Chem. Int. Ed. Engl. 43, 
2117 (2004). 

7 L. Neel, Ann. Geophys. 5, 99 (1949). 

8 W. F. Brown, Phys. Rev. 130, 1677 (1963). 

9 W. T. Coffey, D. S. F. Crothers, J. L. Dormann, Yu. P. 
Kalmykov, , and J. T. Waldron, Phys. Rev. B 52, 15951 
(1995). 



10 L. Neel, C. R. Acad. Science 224, 1550 (1947). 

11 E. C. Stoner and E. P. Wohlfarth, Philos. Trans. London 
Ser. A 240, 599 (1948), reprinted in IEEE Trans. Magn. 
MAG-27, 3475 (1991). 

12 A. Thiaville, J. Magn. Magn. Mat. 182, 5 (1998). 

13 A. Thiaville, Phys. Rev. B 61, 12221 (2000). 

14 W. Wernsdorfer, E. Bonet Orozco, K. Hasselbach, 
A. Benoit B. Barbara, N. Demoncy, A. Loiseau, D. Boivin, 
H. Pascard, and D. Mailly, Phys. Rev. Lett. 78, 1791 
(1997). 

15 J. Kurkijaxvi, Phys. Rev. B 6, 832 (1972). 

16 L. Gunther and B. Barbara, Phys. Rev. B 49, 3926 (1994). 

17 A. Garg, Phys. Rev. B 51, 15592 (1995). 

18 M. Murugesu and et al., to be submitted (2005). 

19 L. Sorace, W. Wernsdorfer, C. Thirion, A.-L. Barra, M. 
Pacchioni, D. Mailly, and B. Barbara, Phys. Rev. B 68, 
220407(R) (2003). 

20 W. Wernsdorfer, N. E. Chakov, and G. Christou, Phys. 
Rev. B 70, 132413 (2004). 



W. Wernsdorfer, M. Murugesu, and G. Christou, 
cond-mat/0508437 . 

W. Wernsdorfer, E. Bonet Orozco, K. Hasselbach, A. 



Benoit, D. Mailly, O. Kubo, H. Nakano, and B. Barbara, 
Phys. Rev. Lett. 79, 4014 (1997). 



